What is Claimed Is: 

1 1. A method for decomposition of projection data acquired by scanning a set of objects 

2 using at least two x-ray spectra, the projection data including low energy projections (P L ) and 

3 high energy projections (P H ), said method comprising: 

4 A. solving the projections P L and P H to determine a photoelectric line integral (A p ) 

5 component of attenuation and a Compton line integral (A c ) component of 

6 attenuation of the set of scanned objects using multi-step fitting; and 

7 B . reconstructing a Compton image 7 C and a photoelectric image 7 P from A c and A p . 

1 2. The method of claim 1, further including, prior to step A, performing a calibration 

2 procedure using simulated data or measured data or some combination of simulated and 

3 measured data. 



13. The method of claim 1, further including, prior to step A, performing a calibration 

2 procedure, wherein said calibration procedure includes generating low energy iso-transmission 

3 contours for known values of P L at known values of A p and at known values of A c . 

1 4. The method of claim 3, wherein said calibration procedure includes, for each of said low 

2 energy iso-transmission contours, fitting A p to a polynomial function gi(A c ), wherein g L is a 

3 polynomial function that represents the shape of the contour. 

1 5. The method of claim 4, wherein g L includes a set of coefficients g Li determined at said 

2 known values of P L and said calibration procedure includes fitting the values of each coefficient 

3 g Li to a polynomial function h u (Pi). 

1 6. The method of claim 1 , wherein said calibration procedure includes computing minimum 

2 and maximum values of P H for each of said low-energy iso-transmission contours as a function 

3 ofP L . 
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1 7. The method of claim 6, wherein said calibration procedure includes fitting the minimum 

2 values of P H to a polynomial function m L (P L ). 

1 8. The method of claim 6, wherein said calibration procedure includes fitting the maximum 

2 values of P H to a polynomial function n L (P L ). 

1 9. The method of claim 1, wherein said calibration procedure includes generating high 

2 energy iso-transmission contours for known values of P H at known values of the A p and at known 

3 values of the A c . 

1 10. The method of claim 9, wherein said calibration procedure includes, for each of said high 

2 energy iso-transmission contours, fitting A p to a polynomial function gu(A c ), wherein g u is a 

3 polynomial function that represents the shape of the contour for a given P H . 

1 11. The method of claim 10, wherein said g H includes a set of coefficient g Hi determined at 

2 said known values of P H and said calibration procedure includes fitting the values of each 

3 coefficient g H i to a polynomial function h Hi (P H ). 
4 

5 12. The method of claim 1, wherein step A includes generating a low energy iso-transmission 

6 contour corresponding to P L and a high energy iso-transmission contour corresponding to P H . 

1 13. The method of claim 12, wherein step A includes determining the values of the A p and A c 

2 at the intersection of said low energy iso-transmission contour and the high energy iso- 

3 transmission contour. 

1 14. The method of claim 13, wherein the intersection of the low energy iso-transmission 

2 contour and the high energy iso-transmission contour is determined by equating a first 

3 polynomial function gi(A c ) representing said low energy iso-transmission contour, wherein g L is 

4 a polynomial function that represents the shape of the contour of P L , with a second polynomial 
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5 function gu(A c ) representing said high energy iso-transmission contour, wherein g H is a 

6 polynomial function that represents the shape of the contour for a given P H . 

1 15. The method of claim 12, further including computing a modified value of the input low 

2 energy projection data (P Lc ) and a modified value of the input high energy projection data (P Hc ), 

3 wherein each of said modified values P Lc and P Hc are clamped to be bounded between two 

4 values. 

1 16. The method of claim 15, including representing the low energy iso-transmission contour 

2 with a polynomial function g L and determining a set of coefficients of g L as a function of Plc- 

1 17. The method of claim 15, including representing the high energy iso-transmission contour 

2 with a polynomial function g H and determining a set of coefficients of g H as a function of P Hc . 

1 18. The method of claim 15, further including, prior to step A, generating calibration data 

2 using P u wherein P Lc is computed by clamping the value of P L to lie between 0 and the 

3 maximum value of P L used to generate said calibration data. 

1 19. The method of claim 15, wherein the modified value P Hc is determined by clamping the 

2 value of P H to lie between a minimum value of P H (P H min) and a maximum value of P H (P H max). 

1 20. The method of claim 19, including determining P H min as a function of P Lc and a 

2 polynomial function n L and determining P H max as a function of P^ and a polynomial function m L , 

3 wherein m L is a polynomial function that determines P H min for a given value of and wherein 

4 « L is a polynomial function that determines P Hm m for a given value of P Lc . 

1 21. The method of claim 1, wherein step A includes calculating a scaled Compton line 

2 integral value (^ cs ) as a function of a scale factor s c and A c and calculating a scaled photoelectric 

3 line integral value (^ ps ) as a function of a scale factor s p and A ps . 
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1 22. The method of claim 21, wherein step B includes constructing said I c and said I p as a 

2 function of said A cs and said A ps . 

1 23. The method of claim 1, further including, after step B, determining an image of a basis 

2 function X(/ x ) and a basis function Y(/ Y ), by solving I c and I p on a pixel-by-pixel basis, wherein 

3 the basis functions X(/ x ) and Y(/ Y ) are functions linearly combined to determine the pixel 

4 intensities in I c and 7 P . 

1 24. A method for decomposition of projection data acquired by scanning a set of objects 

2 using at least two x-ray spectra, said projection data including low energy projection data (P L ) 

3 and high energy projection data (P H ), said method comprising: 

4 A. performing a calibration procedure using at least some simulated data or 

5 measured data or a combination of simulated and measured data, including: 

6 i. generating low energy iso-transmission contours for known values of P L 
1 and high energy iso-transmission contours for known values of P H ; 

8 ii. generating a polynomial g L that represents the shape of the low energy iso- 

9 transmission contour for each P u wherein g L includes a set of coefficients 

1 0 g u determined at said known values of P L ; 

1 1 iii- generating a polynomial g H that represents the shape of the high energy 

12 iso-transmission contour for each P H » wherein g H includes a set of 

13 coefficients g H{ determined at said known values of P H ; 

14 iv. generating polynomials h L that represents the variation of the coefficients 

1 5 of the polynomial g L as a function of P L ; 

16 v. generating polynomials h H that represents the variation of the coefficients 

1 7 of the polynomial g H as a function of P H ; 

18 vi. determining the minimum and maximum values of P u for each 

1 9 transmission line corresponding each P L \ 

20 vii. generating a polynomial m H that represents the variation of the minimum 
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2 1 value of P R as a function of P L \ and 

22 viii - generating a polynomial n H that represents the variation of the maximum 

2 3 value of P H as a function of P L \ 

24 B - solving the projections P L and P H to determine a photoelectric line integral (A p ) 

25 component of attenuation and a Compton line integral (A c ) component of 

26 attenuation of the set of scanned objects using a multi-step fitting procedure, 

27 including: 

28 i. computing the values of each coefficient gu using a polynomial function 

29 h u (PO and computing the values of each coefficient g H i using a 

30 polynomial function AhiCPh); and 

3 1 ii- determining v4 c and ^ p as a function of P L and P H , using the coefficients of 

32 g L and the coefficients of g H ; and 

33 C reconstructing a Compton image I c and a photoelectric image I p from ^ c and A p . 

1 25. The method of claim 24, further including, after step C, determining an image of a basis 

2 function X(/ x ) and a basis function Y(/ Y ), by solving image I c and image 7 P on a pixel-by- 

3 pixel basis. 

1 26. A system for decomposing projection data for a set of scanned objects acquired using at 

2 least two x-ray spectra, said system comprising: 

3 A. media for storing low energy projection data (P L ) and high energy projection data 

4 TO; 

5 B. a decomposition module configured to determine a photoelectric line integral (A p ) 

6 component of attenuation and a Compton line integral (A c ) component of 

7 attenuation for P L and P H using multi-step fitting; and 

8 C. an image construction module configured to construct a Compton image (I c ) and a 

9 photoelectric image (7 P ) from the A p and A c . 



1 



27. The system of claim 26, wherein the decomposition module includes: 
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2 D. a calibration module configured to calibrate the decomposition module using at 

3 least some simulated data or measured data or a combination of simulated data 

4 and measured data. 

1 28. The system of claim 27, wherein the calibration module is configured to generate low 

2 energy iso-transmission contours for known values of P L at known values of A p and at known 

3 values of A c . 

1 29. The system of claim 28, wherein the calibration module is configured, for each of said 

2 low energy iso-transmission contours, to fit A p to a polynomial function gJA c ) 9 wherein g L is a 

3 polynomial function that represents the shape of the contour. 

1 30. The system of claim 29, wherein g L includes a set of coefficients g Li determined at said 

2 known values of P L and the calibration module is configured to fit said set of coefficients g Li to a 

3 polynomial function h u (Pl). 

4 31. The system of claim 28, wherein the calibration module is configured to compute the 

5 minimum and maximum values of P H for each of the low-energy iso-transmission contours 

6 corresponding to P L . 
7 

1 32. The system of claim 28, wherein the calibration module is configured to fit the minimum 

2 values of P H to a polynomial function m L (P L ). 

1 33. The system of claim 28, wherein the calibration module is configured to fit the maximum 

2 values of P H to a polynomial function n L (P L ). 

1 34. The system of claim 27, wherein the calibration module is configured to generate high 

2 energy iso-transmission contours for known values of P H at known values of A p and at known 

3 values of A c . 
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1 35. The system of claim 34, wherein the calibration module is configured, for each of said 

2 high energy iso-transmission contours, to fit A p to a polynomial function of gw{A c ), wherein g H is 

3 a polynomial function that represents the shape of the contour for a given P H . 

1 36. The system of claim 34, wherein g H includes a set of coefficients g H{ determined at said 

2 known values of P H and the calibration module is configured to fit said set of coefficients g H \ to a 

3 polynomial function h H i (Ph). 

1 37. The system of claim 26, wherein the decomposition module is configured to generate a 

2 low energy iso-transmission contour corresponding to P L and a high energy iso-transmission 

3 contour corresponding to P H . 

1 38. The system of claim 37, wherein the decomposition module is configured to determine 

2 the values of A p and A c at the intersection of said low energy iso-transmission contour and the 

3 high energy iso-transmission contour. 

1 39. The method of claim 37, wherein the decomposition module is configured to compute a 

2 modified value of the input low energy projection data (P^) and a modified value of the input 

3 high energy projection data (P Hc ), and configured to clamp said modified values P^ and P Hc 

4 between two values. 

1 40. The system of claim 39, wherein the decomposition module is configured to clamp the 

2 values of P L to lie between 0 and the maximum value of P L used to generate a set of calibration 

3 data and to compute the modified value P^ as a function of the clamped values of P L . 

1 41. The system of claim 39, wherein the decomposition module is configured to clamp P Hc 

2 between a minimum value of P H (Pumin) and a maximum value of P H (^Hmax). 
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42. The system of claim 43, wherein the decomposition module is configured to determine 
^Hmin as a function of P u and a polynomial n L and to determine P H max as a function of Pu and a 
polynomial m u wherein m L is a polynomial function representing the coefficients of P L for the 
minimum values of P H and wherein n L is a polynomial function representing the coefficients of 
Pi for the maximum values of P H . 

43. The system of claim 26, wherein the decomposition module is configured to calculate a 
scaled Compton line integral value (^ cs ) as a function of a scale factor s c and A c and to calculate a 
scaled photoelectric line integral value (A ps ) as a function of a scale factor s p and A ps . 

44. The system of claim 45, wherein the image reconstruction module is configured to 
reconstruct said I c and said 7 P as a function of said A cs and said A ps . 

45. The system of claim 26, wherein the image reconstruction module is configured to 
determine an image of a basis function X(/ x ) and of a basis function Y(/ Y ), by solving I c and / p 
on a pixel-by-pixel basis, wherein the basis functions X(/ x ) and Y(/ Y ) are functions linearly 
combined to determine the pixel intensities in I c and 7 P . 

46. A system for decomposing projection data for a set of scanned objects acquired using at 
least two x-ray spectra, said system comprising: 

A. media for storing low energy projection data (P L ) and high energy projection data 

TO; 

B. a calibration module configured to calibrate the decomposition module using at 
least some simulated data or measured data or a combination of simulated and 
measured data, and configured to: 

i. generate a low energy iso-transmission contour corresponding to P L and a 
high energy iso-transmission contour corresponding to P H ; 

ii. generate a polynomial g L that represents the shape of the low energy iso- 
transmission contour, wherein g L includes a set of coefficients g u 
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1 2 determined at said known values of P L ; 

13 i»- generate a polynomial g H that represents the shape of the high energy iso- 

14 transmission contour, wherein g H includes a set of coefficients g Hi 

1 5 determined at said known values of P H ; 

1 6 iv - generate polynomials h L that represents the variation of the coefficients of 

1 7 the polynomial g L as a function of P L ; 

1 8 v - generate polynomials h H that represents the variation of the coefficients of 

1 9 the polynomial g H as a function of P H ; 

20 vi. determine the minimum and maximum values of P H for each transmission 

2 1 line corresponding each P L ; 

22 vii. generate a polynomial m H that represents the variation of the minimum 

23 value of P H as a function of P L ; and 

24 viii. generate a polynomial « H that represents the variation of the maximum 

25 value of P H as a function of P L ; 

26 C. a decomposition module configured to determine a photoelectric line integral (A p ) 

27 component of attenuation and a Compton line integral (A c ) component of 

28 attenuation for P L and P H using multi-step fitting, and configured to: 

29 i- compute the values of each coefficient g u using a polynomial function 

30 h u {P\) and to compute the values of each coefficient g Hi using a 

3 1 polynomial function hmiPn); and 

32 ii- determine A c and A p as a function of P L and P H using the coefficients of g L 

33 and the coefficients of g H ; and 

34 D - an image reconstruction module configured to reconstruct a Compton image (I c ) 

35 and a photoelectric image (7 P ) from the ^ and A c . 

1 47. The system of claim 46, wherein the image construction module is configured to 

2 determine an image of a basis function X(Y X ) and of a basis function Y(/ Y ), by solving image I c 

3 and image I p on a pixel-by-pixel basis. 
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